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ABSTRACT 

Main results of this paper are the formulas for the mean 
convection velocities, u, of a viscous fluid and for the mean 
temperature difference in the bulk of the convecting fluid. 

These have been obtained by several ways: by scaling analysis 

of the Boussinesq equations, by analysis of the energetics of the 
process, and by using similarity and dimensional arguments. The last 
approach defines the criteria of similarity and allows the 
proposition of some self-similarity hypotheses. By several 
simple new ways an expression for the efficiency coefficient 
Y of the thermal convection in transforming the heat supplied 
to the fluid into the rate of generation of kinetic energy 
of convection has also been obtained. This expression coincides 
with the one obtained earlier by Lloboutry ( 1972 ) and 
Hewitt ^ (1975) , though these authors gave 

different meaning to this coefficient. An analogy is point€*d 
out between non-turbulent convection of a viscous fluid 
and the structure of turbulence for scales less than 
Kolmogorov’s internal viscous microscale of turbulence. The 
formulas obtained for a and y agree quantitatively well with the 
results of numerical calculations of convection in the Earth’s 
upper mantle (McKenzie et al. 3 1974; Hewitt at al. , 1975; Houston 
and De Bretnaeeker, 1975) which x^ere carried out for a broad 
range of values of the heat flux introduced into the mantle. 
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It is found that the formula for the convection mean velocity 
also describes well some results of numerical calculations by 
Williams (1967, 1971) who studied the flow structure in rotating 
annuli, where atmospheric general circulation was modelled. In 
those experiments the Reynolds number was of order 100. In order 
to clarify the limits of applicability of our fomula for u special 
simple experiments have been performed. In these experiments, I 
checked the proposed theoretical linear dependence of u on the 
depth of convecting fluid. The measurements revealed that this 
dependence is observed with a fair accuracy at least up to 
Re-vlO^. 

^ext, a simple phenomenological theory of gravitational 
convection is developed when motions arise due to density 
diff erentation on the lower boundary of fluid or within it. 

It is shown that this case and the thermal case can be 

described in a unified way by introducing the density deficit 
flux. Applications of the theories to mantle 

convection are discussed rather briefly. The observed value - 

of the geothermal flux gives a limitation to the intensity 

of gravitational convection and to the rate of differentiation 

of the mantle material. A simple experiment has been devised and 

its results confirm some basic conclusions of this elementary 

theory. Applications of these theories to mantle convection 

are discussed briefly. 




/ 

In the last section an attempt is made to classify some 
forced flows. The two types of convection described here and 
turbulence beyond the Kolmogorov microscale are one type of 
flow. For these kinds of flows the total kinetic energy does 
not depend on the mass of the fluid in motion. The last 
property also exists in the theory of similarity for the 
circulation of the planetary atmospheres (Golitsyn, 1970, 

1973). The turbulent convection, in the boundary layer, or 
in the fixed volume (described in Appendix 3) , some models 
of general circulation of a non— rotating atmosphere (Leovy 
and Pollack, 1973; Burangulov and Zilitinkevich, 1976) and 
the turbulence in the inertial range of scales form another 
type of flows which could be called inertial or turbulence 
dominated. All of these aforementioned flows (or flow 
models) have one common feature, obeying "the rule of the 
fastest response," which says that the total kinetic energy 
of the fluid system is equal to the total supplied power 
multiplied by the shortest relaxation time of the system. The 
rule, found empirically by the author, is not a universal one 
but it may have some pragmatic value. 

Appendix 1 is devoted to a discussion of the importance of 
the convection efficiency coefficient. In Appendix 2 the 
connection of the elementary theory developed here with the usual 
description of convection in particular in terms of the Rayleigh 
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number, Ra. It is shown that for large Ra the temperature field 
is almost isothermal within the bulk of fluid and the main 
changes are within the thermal boundary layers . 


INTRODUCTION 


Convection of a viscous fluid is discussed in 
thousands of papers, dozens of books and reviews and nevertheless, 
due to the nonlinearity of the phenomenon, the problem is still 
far from a complete solution and each new setting of the problem 
requires a new analytic, numerical or experimental approach. 
Theoretically and experimentally convection has mainly been 
studied by prescribing the temperature difference at the 
surfaces bounding the fluid. Much fewer are cases 
when the heat flux was given >as is the case for most 
natural phenomena where convection is observed. The mere 
change in the setting of the problem allows one to obtain 
by elementary ways a number of sufficiently general and useful 
results described here. 

This paper originated from my old interest in the problem 
of motions in the Earth’s upper mantle and only quite recently 
I realized its more general meaning. Due to existing large 
uncertainties in our knowledge on the mantle properties, in 
the details and even in the causes of the convection in the mantle 
my goal was only to obtain simple estimates of mean velocities in 
dependence upon the determining parameters and by simple but 
general ways . 

First preliminary estimates of velocities within the mantle 
were obtained in 1970 (in an unpublished paper) using some 
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"principles'* noted by me in the theory of similarity for circulation 
of planetary atmospheres (Golitsyn, 1970, 1973), However, due to 
some reasons which became clear to me only recently, those 
estimates gave velocities an order of magnitude larger 
than the observed several centimeters per year — • velocities of 

lithospheric plates. An important paper by McKenzie ^ (1974) 

which I learned of only in the fall of 1976 due to R, Hide, served as 

an impulse to turn back to the Droblem, In the 

paper a detailed account of plate tectonics is presented 

(see also Turcotte, 1975), some estimates of convection 
characteristics are given on the basis of boundary layer 
arguments and, which is most impressive, the results of numerous 
numerical experiments are described in modelling convection in 
the upper mantle. Impressed by this paper I developed a 
similarity theory for the convection in a very viscous fluid 
which gave by dimensional arguments, formulas which agreed well 
with the numerical results by McKenzie ^ (1974) * This 

stage of work was described in my short note (see Golitsyn, 1977a). 
While this note was in press, I wrote a letter to D. P. McKenzie 
who kindly sent me reprints of many of his papers. The most 
useful for me was the paper by Hewitt, McKenzie and Weiss (1975) 
on dissipative heating in convective flows. The results of the 
paper allowed me to <^heck some conclusions of my first note and 
presented the possibility to understand the energetics of many other 
types of convective flows, in particular, of experiments in rotating 
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annuli on modelling atmospheric general circulation (this 
particular problem is considered in Golitsyn, 1977b). All this 
has allowed me to develop a rather general approach to the 
problem of viscous fluid convection which has a broader 
applicability than to flows with small Reynolds number and very 
large Prandtl number as it had been thought at first (Golitsyn, 
1977 a). As a by-product I noted some general features and 
similarities for a number of forced flows. 


The present paper is constructed in the following way. 
First by the scaling analysis of the convection equations the 
estimate of the mean velocity: 



is obtained for a cell with the size d. Here a is the thermal 
expansion coefficient, g is the gravity acceleration, f is the 
density of the heat flux > Ufs the dynamic viscosity, and is 
the specific heat at constant pressure of the fluid. The 
analysis also gives a mean temperature difference in the buld 
of the convecting fluid: 



where v -y/p is the kinematic viscosity and p is the density of the 
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fluid- Then an estimate is given of the efficiency coefficient y 
of convection in transforming the supplied heat into the rate of 
generation of kinetic energy which — the rate — in a 
steady state is equilibriated by the mean rate of dissipation. 
The estimate is 

d agd 

^ " T = [3] 

These results ware evidently first obtained by Lliboutry 
(1972) from different arguments and also by Hewitt £t (1975) 
by analysis of energy and entropy balance (see also 

Golitsyn, 1977b). The other derivation of the formula for u 
if Y is -known may be obtained from the expression for 

the rate of kinetic energy dissipation per unit mass of a 
viscous fluid. This derivation specifies the value of the 
numerical coefficient in the formula [1] and [2]. 

Then an analysis of the similarity criteria is carried out 
for the convection of a viscous fluid. It is found that if we 
scale velocities by [1] and temperature by [2] then the equations 
do not contain any numerical values for small Reynolds numbers 
(except in the boundary conditions) which means that such 
convective flows are self-similar at least in tie bulk of fluid 
(if the one length scale is present) . Another type of dimensional 
argument is also produces the formula [1]. 


An exact analogy is noted between the convection of a viscous 
fluid and the structure of a velocity field in the developed 
turbulence for scales less than the Kolmogorov viscous 
microscale. To make the analogy more convincing the Kolmogorov 
formulas are obtained also by dimensional arguments which require 
some slight reformulation of the problem. Then it is shown that 
the analogy is valid for any type of non-turbulent convection 
disregarding the value of the Prandtl number. 

Then the results are compared with the data from 

various numerical experiments un modelling the mantle convection 

(McKenzie et al. 1974; Hewitt 1975; Houston and De Bremaecker, 

1975), and on modelling convection in annuli (Williams,- 1967, 1971). 

In order to understand the applicability of the formula [1] to the 

convection with a rather high Reynolds number, a few simple 

experiments were carried out which show that the expression [1] is 

q 

valid up to Re^lO at least. Some discussion of this fact is 
presented. 

Then gravitational convection is considered X'^hich takes 
place due to density differentiation of the material at the 
boundary or within the fluid. Such convection in many respects 
is found to be similar to the thermal case. Moreover, it is possible 
to describe both types of convection in a uniform way by introducing 
the flux of the density deficit. Simple experiments confirm some 
conclusions of the theory. 
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After that there are applications of the results to the 
Earth's upper mantle convection and questions of its modelling 
are discussed. In the last section an attempt to classify some 
geophysical forced flows is made and a "rule of the fastest 
response" is formulated. 

Though some of the results (e.g., [3]) are known or 
practically known (e.g,, [1] may be readily obtained from the 
formulas [47] and [53] of McKenzie , 1974) , they are 

obtained here by another and direct way which clarifies their 
meaning and shows the limits of their applicability. Other 
theoretical and experimental results are new. 

The discussion of all of this material in a single paper 
seemed to me useful for the construction of a uniform view 
various phenomena quite different at first sight. 

2 . Scaling analysis and energetics of convection . 

Convection of an Incompressible fluid is usually described 
by the Boussinesq equations : 

Bt; n 

1 0 

TTT )z; = Vp + vAi? [4] 

9 k k. P ^ 

V- y ~ 0 


[5] 
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In the momentum equation [4] t? is velocity, p is pressure, 

0 

p is density, T is the temperature deviation from its equilibrium 
(adiabatic) distribution in the absence of convection. In the 
energy equation [6] q is the intensity of heat sources per unit 
mass and £ is the heat intensity due to friction. Boundary 
conditions must be added here of no-slip for velocities, the heat 
flux on the lower boundary and say, the temperature at the upper 
level (for the problem of mantle convection) . 


An analysis of this set of equations and of their applicability 
to the upper mantle convection has been given by McKenzie ^ a^. (1974) . 
When the mantle is further considered I shall use their values 

3 3 3 

of the parameters: p“= 3,7 • 10 kg/m , = 1.2 • 10 J/kgK, 

2 -5—1 17 2 

g = 10 m /sec , a = 2 *10 K , V= 2 *10 m /sec, 

— 6 2 

k - 1. 5 • 10 Tti /sec. The mean value of the geothermal flux f 

2 2 — 6 2 
is about 0.06 W/m = 60 erg/cm sec = 1.4 • 10 cal/cm sec. For 

the mantle conditions the left-hand side of equation [4] is absolutely 

insignificant because the Reynolds number for a layer depth of 

5 —Q 

d = 7 • 10 m = 700 km and velocities of order 10 m/sec = 3 c»n/yr 
'-20 

is of order 10 
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For further analysis it is convenient to use instead of [4] the 
equation for vorticity w = curl V and instead of temperature the 
entalpy difference e]^ - e - e^ = Cp(T - Tq) , where Tq is the tempera- 
ture of the upper boundary. Then instead of [4] and [6] we shall 
have for Re « 1 

vAw = H X n, n = g/g [8] 

8e]^/3t + yVe^^ = kAej^ + g + e, [9] 

where H = c^/ag [10] 

is a quantity with the dimension of length. If the medium is a 
perfect gas then g/cp = (dT/dz)^^^ is the adiabatic lapse rate, a = 1/T 
and H = T/(dx/dz) and i.e., it is a depth of adiabatically stratified 
gas layer (for air H - 30 km) . For other media H may be considered 
as a characteristic depth of the layer of material stratified by 
gravity field. For the upper mantle H = 6000 km. For water 
a = 2»10"'^K"7^ Cp = 4.2*103j/kgK and H = 2000 km. 

Further on for simplicity we shall consider convection with 
large Peclet numbers (though this requirement is not quite necessary -- 
see the end of §4) when the diffusional heat transfer is small as 
compared with the convection one. For this the Prandtl number of the 
fluid Pr = V/k = Pe/Re should be large enough. For the mantle»e.g.» 
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Pe 500 because Pr lO^^. As a rule the contribution of 
viscous heating in eq. [9] is also small (conditions for its neglect 
will be discussed later; the special study of its role in this 
respect have been carried out by Hewitt et al. , 1975). As a result 
in the stationary case the eq. [9] becomes 

yVei q . [9'] 

The total heat flux observed at the upper boundary is equal to 

/ = pqd + . [11] 

Here for simplicity after McKenzie et al. (1974) we consider 
q= const, otherwise the difference betvjeen fluxes at two boundaries 
of the layer is / *^1 “ It is convenient to introduce 

again after McKenzie et al. (1974) 

3 = pqd// = pqd/(pqd + [12] 

the part of the heat flux produced within the layer to the total flux. 

For 3 = 0 the heat is introduced only at the lower boundary and for V 

3 = 1 the heat is produced only within the layer. The value of q 

could be estimated as | 

q = 3//pd a; //pd [13] I 


i 
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The simplest estimate of the characteristic velocity of con- 
vective motions U can be obtained from the scale analysis of eqs, [8] 
and [9']. For the case of a bounded cavity with one characteristic 
length scale d (note that computations by McKenzie et al- , 1974 and 
Hewitt et al. , 1975 were done for the square region) we have from 

[9*] accounting for [13] that Ve^^ ^ q/U ^ //Upd. Substituting this 

2 q 

into [8], taking into account that 0)/d U/d*^ we may finally 

get 



For the characteristic temperature difference we get from [9’] 

then 



Because [9’] disregards thermal boundary layers (term kAT is 
omitted) this value is evidently scaling the temperature difference 
within the bulk of fluid, i.e., outside the boundary layers. For 
special discussion of the temperature field in the convective fluid 
see Appendix 2. 

From eq. (7) for the rate of energy dissipation e \)uVd^. 
Substituting here U from [1*] we get 

e ^ qd/H (//pd)(d/H), [14] 
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The quantity 


Y = d/H = agd/Cp [3] 

can be named the efficiency coefficient of convection in transforming 
the supplied heat into the rate of generation of kinetic energy of 
convective motions. This rate of generation in the steady state must 
la equaled in average by the rate of the kinetic energy viscous dis- 
sipation. 

We may neglect e in the energy equation [9] 
if Y for the upper mantle y 700 l<m/6000 km = 0.117. It is 

because of this low value that the computations by McKenzie and his 
colleagues with and without viscous heating produced practically 
the same results on the patterns and intensity of convection. Hewitt 
et al. (1975) have shown that integrally, as averaged in space and 
time the viscous dissipation is exactly equal to the work of pressure 
forces and this fact clarifies why viscous dissipation is equal in 
average to the rate of generation of kinetic energy of convection. 

In the same paper by Hewitt et al. (1975) there is a derivation 
of the formula [3] for the efficiency Y based on consideration of 
the global balance of energy and entropy in the fluid. With some 
corrections insignificant for the final result such a derivation is 
reproduced also by Golitsyn (1977b) . However, the formula [3] was 
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evidently first obtained by Lliboutry (1972). He used very 
simple arguments deserved for reproduction here. When a unit volume, 
which is lighter by p' = paT’ of its environment, is lifted on the 
height d then a potential energy paT*gd is gained which is transformed 
into heat. An amount of heat carried by the volume during 
convection is CppT* if the diffusion of heat is small. The efficiency 
of convection as a heat engine may be defined as a ration of these 
two quantities 

_ pgT’gd qgd ^ 

^ CppT' - Cp H • 

This derivation is essentially correct though not all the neces- 
sary details have been Included, because convection is not only a heat 
engine but also a motion engine. It shows that the proportionality 
coefficient between e and //pd is exactly unity (if all the heat comes 

from below; if 3 < 1 then Hewitt et al. , 1975, have shown that 
Y = 1 - 3/2 for g = const). The other details are in a derivation 
similar to this one after our formula [41]. At the end of Section 
4 it will be shown that the requirement Pe >> 1 is not essential and 
the estimate y d/H as well as eq. [1] for the mean velocity are 
valid at any Prandtl number if Re is not too large. <?ee also Appendix 3.) 

If we take the expression for the efficiency y as a known, e,g., 
from the above described considerations, then the formula for the 
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mean convection velocity can be determined with accuracy up to a 
multiplier of order unity by a second means. Suppose than in the 
expression [7] for the mean rate of energy dissipation e all 
derivatives are of the same order. This is supported by inspection 
of graphs of the stream function in computations by McKenzie et al. 
(1974). The graphs have a shape of more or less concentric circles 
spaced about uniformly. Approximate as 2U/d. For the plane 

case as in McKenzie et al. (1974) or for roll convection in the 
sum [7] there are eight terms, therefore e ^ 32vU^/d^. At the other 
hand due to [14] e Y//pd ~ //pH (note a rather unexpected fact 
that the dissipation rate is determined only by the heat flux and 
the fluid properties, but not by geometry). * , Equate these 

two expressions for e and get for U: 




/ agf V V 

\32\mJ \32\ic^j 


d 


[ 1 '] 


If the motion is essentially three-dimensional than in the sum 
[7] there are 18 terms and for the same approximation of derivatives 
one has in the eqn. [I’j 72 instead of 32. Therefore all conditions 
being the same the intensity (velocity) of the 3D-convection is not 
more than (72/32)^ =1.5 times less than for 2D-case. 


This derivation of the mean convective velocity is quite general. 
It requires only that viscosity be controlling the flow structure, 
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i.e., the Reynolds number should not be too small and the con- 
vection should not be fully turbulent. The discussion of this 
statement will be continued afthe end of Section 5. 

3. Similarity theory for the convection of viscous fluid . 

The similarity theory is a base for modelling of various phenomena. 
In several cases the formulation of some self-similarity hypotheses 
gives a possibility to get by dimensional arguments some non-trivial 
results. In our case all this allows one to derive sufficiently 
rigorously eq. [1] again and by this to determine the conditions 
for its validity and for modelling convection. 

Let us formulate our problem again for variables w - V x y and 
- Cp(T - Tq) for a general case. Then we have 

dw 

^ - M)v = - I [Vej^ X n] + vA“ [15] 

— = kAej^ + q + e [16] 

with the zero velocity for all boundaries and e^^ s o at z - d and 

kDe^^/Sz = at z = 0 . [17] 


I 


1 
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Transform this system to a non-dimensional form by taking d 
as a length scale, de = (/vH/pd^)"'® from [2] as the enthalpy scale, 
eq. [1] as the velocity scale. Then the time scale is defined as 


T = d/U = (yc„/ctg/)^ = (pH//)'"' . [18] 


It is interesting that the characteristic time of convection, 
i.e. , the fluid particle turnover time does not depend on the geometric 
scales but only on heat flux and fluid properties. The equations 
with all the dependent and independent variables being non-dimensional 


Re[doydt - (Ol)V)i?] = n] 4- AO) 


[15'] 


de;j^/dt == Pe ^Ae^ + q 4- 


[16'] 


3e-|/8z = -Pe at z = 0 and e^ - 0 at z = 1 [17’] 


Here Re = Ud/v is the Reynolds number. It can be written noting 


[18] as 


Re = d2/vT = T /t 


i.^e, , to represent as a ratio of the two relaxation times; of the 
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viscous time = d^/v to the turnover time T = d/U. In the 
energy equation [16’] there are tx^^o similar criteria: first is the 
F4clet number Pe == Ud/k and the second y = cL/H is the convection 
efficiency. 

An important class of convective motions is v/hen all the 
sitnilarity criteria are small or large compared to unity. For the 
upper mantle, e.g., Re f^ 3*10"21j 'Pe ^ 10^, y 10“^. For Pe >> 1 
the relative width of the thermal boundary layer is small and if we 
disregard its structure then we may assume ' self-similarity of 
convection in the main volume of the fluid, because eqs. [15 ’] and 
[16’] will contain none of the similarity criteria, i.e., all the 
terms left will be of the same order. This means that the non- 
dimensional scales of velocity and enthalpy fluctuation are of order 
unity, i.e., estimates [1] and [2] are strictly valid if Re << 1, 

Pe >> 1 and y << 1. At these conditions ^ moreover, one may expect 
the similarity of flow and temperature patterns of convection out- 
side the thermal boundary layer. We shall see later in Sections 4 
and 5 that for validity of [1] and [2] the first two conditions 
are not necessary, the estimates are valid in a much broader range 
of Reynolds and Pdclet numbers though outside the above range one 
should not expect the similarity of flow patterns. 

If beside the vertic.^1 scale d there is a horizontal scale L 
then all the values have to be dependent on the ratio d/L - H, 
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For similarity of convective flow patterns one needs the same 

values of ^ and self-similarity with respect to other criteria. 
The connections of these similarity criteria with Rayleigh numbers 
are discussed in Appendix 2. 

The derivation of eqn. [1] using more traditional similarity 
and dimensional arguments may be found in Golitsyn (1977a). Note 
that there, due to different scaling of velocity and temperature, 
there is another similarity criterium in the momentum equation. 

Still another derivation of eqn. [1] can be proposed if we 
note that our equation system for the case Re << 1, y << 1 and 
q == 0 is invariant relative to a group of transformations (the 
idea of this derivation was suggested to me by G. I. Barenblatt) . 

For this case eqns. [15]-[17] are simplified to 

VHAo) = Vej^ X n [15”] 

de^/dt = kAe^ [16”] 

kde^/^z = ~//p at z = 0 and e^^ = 0 at z = d. [17”] 

The structure of the system allows one to consider the dimension 
of enthalpy as arbitrary, but then the dimensions of external 


' i I I J I I I 
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parameters will be as follows: [b] = [vH] = [e]LT, [f/p] = [e]LT~^, 

[d] = L and [1^1 = L^T"^. Here [e] is the dimension of enthalpy, 

L and T are dimensions of length and time. Note that in eqn. fl5"] 
there is the combination vH = b only and because we neglect inertial 
terms and viscous heating then the kinematic viscosity v separately 
does not enter. 

From the three first external parameters one may construct the 
quantity with the dimension of velocity 

V = (/d2/pvH)^ = (ag//ycp)^d = (f/yH)^d . [1] 


If we account for the fourth parameter, the thermodiffusivity, 
then we may construct from all the parameters one non-dimensional 
combination, which is Pe = Vd/k = (//pH)^(d^/k) . If Pe >> 1 then 
the concrete value of k is not essential and we may not account for 
it in the set of external determining parameters. 

The possibility of an arbitrary choice of the enthalpy 
dimension can be obtained from the next observation (by G. I. Barenblatt). 
The system [15"]-[i7‘'1 does not change its form under transformation 


ej^ ae, vH avH, f/p af/p 


[19] 
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where a is an arbitrary number. If we would consider as usual 

[e] = then [f/p] = L^T"^, [vH] = l3t- 1, [d] = L, [k] ^ l2t-1 

and from the three first parameters one may already construct a 

new non-dimensional criterium II = (//p)^(vH)~3-/^d . Under transforma- 

1/3 

tion [19] this criterium is changing as n a II, i.e. , it 
is dependent on the choice of the value of a, in contradiction 
to the basic system [15"]-[17"j. Therefore, in our case we may 
consider the dimension of enthalpy as arbitrary with total justi- 
fication. 

For a volume with the characteristic size d one obtains the 
formula for the total kinetic energy of the convection using [I']: 

E = i pd%2 = c2 [20] 

2 ^ 2 vH 

where coefficient 1/32 or 1/72 and may also depend on the 

shape of the region. 

The formula [3] for the erlficienGy coefficient y d/H may 
be understood in the following way at least for the case of Re << 1. 

Including dissipation means the appearance of a new dimen- 

sional parameter in the energy equation [16'*], i . e. , the kinematic 
viscosity V itself. Adding it to the existing parameters b = vH, 
f/p and d gives a possibility to produce the non-dimensional similarity 


criterion ~ vd/b - d/H. With respect to dissipation we 
have a typical example of self-similarity of the second kind 
in the terminology of G. I. Barenblatt (1976). The mean convection 
velocity for a viscous fluid is self-similar with respect to II , 

Re and the last independence means that U does not depend on 

V separately, only on the combination b = vH. If we take an 

analogous hypothesis for the dissipation rate e then we could construct 

-1 

the quantity with the dimension of e([e] = [e]T from [16]) from the 
parameters b = vH, //p and d only as e = c^//pd, where would be 
some numerical constant. The neglect of viscous heating means the 
neglect of the similarity criterion y^^ = vd/b. But this can be done 
in determining the dissipation from the dimensional arguments only 
if C]^ = ^l('Yi) const at y^ 0. Here we have, after Hewitt et 
al. (1975) 

llm = Y3_(1 -|) , [21] 

Y1 0 

where g is from [12]. The non-existence of a finite limit for a 
constant when some similarity criterium tends to zero (or infinity) 
is the property of the self- similarity of the second kind. The 
constant usually has a power dependence on the criterium. The 
character of this dependence can not be defined from dimensional 
arguments and some other arguments should be applied (here the 
expression [7] for e has been used) . 


Before seeing the paper by Hewitt et al. (1975) I determined 
the value of y = eM// for the upper mantle empirically (not sus- 
pecting the complicated situation just described) using the results 
of computations by McKenzie et al. (1974). In Figs. 18-20 of that 
paper there are graphs of the horizontal velocity component u at 
the free upper surface of the computational cell in dependence on 
the horizontal coordinate x and the stream functions for three 
types of the lower boundary conditions: 3 = 0, % and 1. Assuming 
as above all derivatives in the sum [7] of the same order (no 
boundary layers for u) we may write e pv(9u/9a:)^. Differentiating 
graphically u(a:) and taking squares of the derivatives we may 
estimate the value of e . For 3=0 (all heat is from below) the 
calculations give (9u/9a:)2 = 0.16 mm^/km^yr^ = 1. 6*10“^^sec~^ because 
1 km = 10^ mm and 1 year = 3.15*10^ sec. From here e = 2.5*10"^^ 
m^sec^^H(W/kg).It needs to be compared with the mean geothermal flux 
/ = 6 10~2wm”2. Then we get y = e// = sM/f - 0.11. The calcula- 
tions by the formula [3] y = d/H give y = 0.117. Hewitt et al. (1975) 
performed special computations of the value of y using their computed 
velocity fields. These computations were done for the mean values 
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of the fluxes through the upper boundary equal to 10 , 10 and 

5.85 Wm“^. In their computations the value of y for 3 = 0 was chang- 
ing from 0.10 up to 0.115. This confirms that the convection 
efficiency does actually depend only very slightly on the heat flux. 
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For g = h the analogous calculations using graphs of Fig. 19 
by McKenzie et al. (1974) give y = 8% and for g = 1 (Fig. 20 ’ — all 
heat from within) I got y 6%. According to [21] we get y = 8.7% 
and 5.8% correspondingly. The exact computations by Hewitt et al. 
(1975) of integral dissipation by velocity fields give changes 
of y for these two cases from 7.3 up to 8.4% and from 4.6 up to 5.4%, 
correspondingly. As we see, the results of the two approximate 
methods and of one computational method of determining the convection 
efficiency are very close between each other. 


Such a value of the convection efficiency coefficient for the 
mantle appeared at first rather high for me* For comparison, the 
similar efficiency for the Earth’s atmosphere in transforming 
solar heat into wind kinetic energy is of order 1%, and for the 
martian atmosphere it is still an order of magnitude less (Golitsyn, 
1973). 


The numerical experiments by McKenzie et al. (1974) 

■ ■ ■■ ■ . . " ■' -4 

were carried out with variation of the heat flux value from 10 up 

to This gives the possibility of checking the dependence 

of velocity on /in the three order range of changes of /, The 
formulas [56], [61] and [62] by McKenzie et al. (1974) describe 
dependencies of computed maximal values of the horizontal velocity 
d on the total heat flux at the upper boundary. 


■ 1 
i 

■ I 

[ 22 ] 


Here Ig = in all Russian literature) , u is measured 

in mm/yr, [/] = Wm*~^. Assuming that the maximal velocity u depends 
on / as the rms velocity u (given by [1]) we see that the character 
of the power dependence of velocity on / as produced by various 
arguments. is fulfilled practically exactly in the numerical experiments 
and does not depend on the way of the introduction of heat into a 
system in accordance with our derivations (say in Section 2). 

If we set in eqn, [20] c = [ (1 "" 3/2)/32]^ as it follows from 
[21] then taking the logarithm of [20] and the numerical values of 
the parameters entering there we may get / in the same units that 

Igu =0.5 Igf +1.76 + 0.5 lg(l - 3/2) ; [23] 

where Ig = 

Remembering that eqns. [22] describe the maximal convection 
velocity, and [23] relates to the rms velocities we see that our 
result is valid not only from the point of view of dependencies on 

Vs 

the external parameters of the problem but it is quite good in deter- j 

mining the value of the numerical coefficient in the formula for 

the mean velocity. | 
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3 = 0 Igu = 0.50 Igf + 1.91, 

3 = % IgQ = 0.49 Igf + 1.80, 

3 = 1 Igd =0.54 Igf + 1.85 . 
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The dependence similar to [22] may be also obtained using the 
results of numerical computations by Houston and DeBremaecker (1975) 
for their case of constant viscosity v = 5* 10^^m^sec“l, density 
p = 3.5*10^ kg/m"“^ and / = 5.94*10“2Wm“2, assuming again that the 
maximal velocity is determined by a formula of the type of [!']• 

The computations give u =16 mm/yr at the distribution of heat 
sources which corresponds to 3 =0.79. The substitution in [1’] of 
the values just given produces 

Igu = 0.5 Igf + 1.82 [22’] 

which agree quite well with eqn. [22]. Some further discussion of 
the results by McKenzie et al. is in Appendix 2. 


4. One useful analogy to convection of viscous fluid 

To clarify the character of mean convection velocity dependence 
on the external parameters let us consider in a similar way one 
classical problem of hydrodynamics but with a slight change of its 
posting to make it and the convection problem alike. Let us con- 
sider the sufficiently developed locally homogeneous and isotropic 
turbulent fluid and a volume in It with a characteristic size d 
smaller than the Kolmogorov viscous microscale 


n = v3/4 ^”1/^ 


[24] 
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For distances r less than ri Kolmogorov (1941) determined the 
structure tensor of the velocity field fluctuation. He found that 

D^^/r) = 2D^j^rr; = 2 r2/15 , [25] 

where D and D are the longitudinal and lateral structure functions, 

II 

the mean squares of the differences of the corresponding velocity com- 
ponents taken at two points separated by the distance r. Kolmogorov 
obtained formulas for D and from the expression for the mean 

It 

rate of kinetic energy dissipation (see [7]) taking into account the 
isotropy and smoothness of velocity gradients at distances r < ri. 

Now we shall obtain an expression analogous to [25] by another 
way not expanding the velocity field in Taylor series on r as 
Kolmogorov did. Let us consider what determines the kinetic energy 
E of the fluid volume with the size d relative to its nearest environ- 
ment at distances also of order d. That means we shall compare the 
kinetic energies of two neighboring volumes with the sizes d < n* 

The relative kinetic energy of such a volume should depend on its 
size d, on the viscosity v and on the total rate of the energy in- 
flow to this volume. 


Q = eM = epd^. 


[Q] = . 


[26] 



Using d, V, and Q we can by this unique way construct the 


quantity with the dimension of the energy 

E = C2Q*d2/v = C2 Qt^ , 


[ 27 ] 


where C 2 is some numerical constant. From here and [26] the value 
of the velocity square in this volume is 



2g2Q 

vpd 



[28] 


The identity of the expressions [25] and [28] is obvious but 
the last one is obtained by dimensional arguments similar to the ones 
used in the preceding section. The structure of these . formulas is 
quite analogous to the structure of the formula [IM for the mean 
velocity of convection of viscous fluid in a bounded volume with 
the size d if one remembers that in the convection e ^ //pH (see 
e.g., [14]). Then the square of the velocity scale from [1] can be 
written as 


.9 e .2 
U V d 


Up to now in order to obtain formulas for the mean velocity and 
efficiency of convection we have assumed Re << 1 and Pe » 1, so 
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Pr > lAi/Rii >> 1, The two last conditions are not essential. 
After the just found analogy it is clear that any laminar convection 
in the sense of its intensity should be similar to the turbulent 
microstructure. In order to show this formally, write the energy 
equation [16] as 


Be 

3t 


+ 




+ e 


Here 



A. 

pCp 


_ k 3e 

pCp 9x^ 



is the kinematic heat flux, where for simplicity the distribution 
density of the volume heat sources q is supposed to be homogeneous 
in the horizontal. Assuming e to be small and the scale of the 
heat flux known independently of its origin we get from here and from 
[15*] the same estimates [1] and [2] for the mean velocity u and the 
temperature difference* 6T. 


’^For the purely diffusive flux /* = k 9T/8z -k6T/d from [2] we may 


get quite naturally 
number. 



kv 

agd^ 


Ra , where Ra is the usual Rayleigh 
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Second, as it follows from the results of experiments described 
in the next section, the condition Re >> 1 is also not necessary 
for the validity of the estimate [1]* It appears that the Reynolds 
number related in the convective conditions to the Rayleigh and 
Prandtl numbers (see Appendix 2) should not be too large in order 
that the convection become not fully turbulent (see Appendix 3). 

Up to the appearance of this regime the analogy between intensities 
of the turbulence viscous microstructure and convection of a viscous 
fluid holds up. 

5. Various experiments on modelling the convection in the 
upper mantle and in other geophysical phenomena 

It has already been noted that due to self-similarity of the 
convection^ for the modelling of the mantle convection it is necessary 
to have Re << 1, and Pe >> 1, The condition y i is always held in 
the laboratory. It is also necessary to have a fluid with Pr =\)/k 
>> 1 large enough. As is shown in Appendix 2 , the Reynolds and 
Pdclet number are related to the Rayleigh number and for setting 
up the convection the last number must exceed its critical value. 
Therefore, in a model it is not necessary to try to reach exact 
coincidence of the corresponding similarity criteria for the model 
and the mantle though the very fact of the self-similarity and the 
limits of its existence should be specially checked in the 
laboratory. 
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Let us consider from this point of view the laboratory experi- 
ments by Booker (1976). For the mantle convection an important 
factor is the dependence of viscosity on depth connected mainly 
with the dependence of viscosity on temperature (see, e.g., Carter, 

1976). Therefore Booker chose as the working liquid a special kind 
of oil with the dynamic viscosity changing from 200 P (poise) at 
-20°C up to 0.2 P at -i-80°C. 

The Rayleigh number Booker defined for the mean temperature of 
the liquid T = %(Ti + T 2 ) where T^^ and T 2 are temperatures of the 
lower warm and of the upper cold surfaces. The structure of the 
convective roll cells depended little on the variation of the liquid 
parameters with height. Using temperature dependencies for the 
parameters written in the Booker paper, we may find for +30°C that 
V = 2. 2*10’’^m^sec*“^, a = 5.7*10'’^K.*“^, k = 7*10”'®m2sec’*^, Cp = 

2.1*10^ /kgK, p = 846 kg/m^. The height of the cell was in the 

limits 1.43 - 1.50 cm. Then the Rayleigh number Ra = agATd^/kv ^ 10^. 
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The Nusselt number in his measurements was found to be Nu = 0.184 Ra * , 

wherefrom Nu ^ 4.6. 

Now by formula 

/ = kpCpATd"^Nu [30] 

we may calculate the heat flux into fluid. For AT - 100 K we get 
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/ = 0.4 Wcm'"^. Having this value we determine the mean convection 
velocity from [1* ] wherefrom V 0.6 mm/sec. Though Booker does 
not report the convective velocity estimates this value seems reason- 
able (from my limited experience of observing convective flows). 

For his experiments now we find Re = 0.04, Pe = 120, and Pr = 3000, 

—8 

and Y = 4*10 

Therefore these experiments satisfy all the conditions of the 
self-similarity for convection. The very large difference in the 
value of the similarity criterion for the dissipation (generation) 
of the kinetic energy = d/H is also not important. The experience 
of the laboratory and numerical modelling of the atmospheric general 
circulation in annuli (see, e. g.. Hide and ilason, 1975, Dolzhansky 
and Golitsyn, 1977, Williams, 1967, 1971, Golitsyn 1977b) and some 
other examples show that the flow pattern does not much depend on 
the value of this criterion. It determines only the intensity of 
the flow. Therefore the experiments of the Booker type allow in 
principle the extraction of much qualitative and even quantitative 
information of the flows in the Earth’s upper mantle. The first 
goal in such experiments should be the check of the dependence 
V ^ or of the whole dependence [1] and [3] . 

Quite surprising was the fact that the dependence of the mean 
convection velocity on the external parameters similar to [1] is 
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also valid for motions in bounded vessels when the Reynolds number 
composed by their rms velocity and the characteristic vessel size 
is not small but has the value of order several tens or even hundreds. 
This had been found while analyzing the energetics of convection of 
water in rotating annuli using the data of detailed computations of 
the process by Williams (1967, 1971). Though the heating and cooling 
from lateral boundaries and rotation decrease the efficiency 
several times in comparison with the heating from below of the non- 
rotating fluid (the situation has been analysed in some detail by 
Golitsyn, 1977b), nevertheless if we know the specific dissipation e 
the mean velocity can be well estimated by the formula 

U (e/32v)^ d [31] 

for the axisymmetric case of convection (Williams, 1967) and by 

u % (e/72v)^ d [32] 

for the three-dimensional case when baroclinic waves are developing 
in an annulus (Williams , 1971). Note the similarity of these formulas 
with Kolmogorov's formulas [ 25) or [28]. 

In fact, using the results of the direct computations by 
Williams (1967) one can find for his case A3 the rms velocity 
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V = 2.6 Tnm/sec. The dissipation e can also be found from his results 
e = 3*10“^ cm^/sec“^. Using now the depth of his vessel d == 5 cm 
one obtains using [31] V = 2.4 mm/sec. For the baroclinic case the 
analogous estimate by [32] is 1.1 mm/sec compared with the computed 
value of 1.2 mm/sec (e = 1.1*10 cm^sec-^, d = 3 cm; for details see 
Golitsyn, 1977b). 

The Reynolds numbers are for the two cases of order 100 and 40, 
correspondingly. One may say that these values are on the one hand 
sufficiently large in order that boundary layers are relatively 
thin but on the other hand the values are still small in the sense 
that the flow pattern is still laminar though rather complicated, 
i.e. , the viscosity forces determine essentially the flow patterns 
and the velocity gradients are appreciable even in the bulk of the 
fluid. 

In an attempt to understand why this theory works even at rather 
high Reynolds numbers, one may formally introduce the Kolmogorov 
internal microscale defined by [24], if the value of e is known. 

For the first William’s case one gets ri = 1.8 mm and for the second 
x\ = 1.3 mm. The ratio of the scale d to r\ is for both cases close 
to 25. Apparently if the value of the criterion 

% = d/n (Re)'2 


[33] 
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is not too large the flow has a laminar or slightly Irregular 
character and our formula [1] may be extended to these conditions,* 
though one must not expect the similarity of the flow patterns for 
Re > 1. The second derivation of our formula for the mean convective 
velocity where the efficiency y and the expression for e are used, 
is also supporting these arguments because it again requires only 
the importance of viscosity for shaping flow patterns and not Re<<l. 

However, Williams has published the detailed results only for 
the two cases. In order to see whether the agreement is not fortui- 
tous I began to think on the possibility of the experimental check 
of eqn. [1] for enlarged Reynolds numbers. I was not able to find 
the necessary data in the literature. The discussion of possibili- 
ties to undertake a check of the full or partial dependencies in 
[1] revealed to me that it is quite a job requiring thorough prepara- 
tion, careful measurements and extensive treatment of data. 

On one of those days (mid-April, 1977) I was asked at home to look 

* Similar results have been obtained by Golitsyn and Steklov (1977) 
while determining the height of the turbopause of a planetary atmo- 
sphere. For the case of Earth the height of the turbopause, i.e. , 
the level where eddy and molecular diffusivities become equal, is 
determined as the height where the microscale rj reaches about 0.1 
of the atmospheric scale height. This was found to be reasonably 
well for other planets too. 
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for the soup boiling on the gas stove. V/atching for grains and 
bubbles in the soup, I estimated that their velocities were of the 
order of several centimeters per second, guessed that the power to 
the kettle 20 cm in diameter was of the order of several hundreds 
of watts and looked what formula [1*] was producing. It gave about 
5 cm/sec. At that moment, I decided to do the necessary experiments 
myself right at the kitchen. My hopes were based on the observations 
of the particles* trajectories. At temperatures before boiling 
the distinct and rather long parts of the trajectories were almost 
rectilinear and horizontal. In space and time the parts (and 
velocities) were rather irregular, which required sufficient 
statistics. 

I realized that I might perform an easy but quantitative check 
of the linear dependence of the mean velocity on the fluid height 
d, keeping other parameters constant, i.e., keeping the fire and 
temperature of the water constant. The experiments took two 
evenings. At first I had worked out the technology and got pre- 
liminary results shox>7ing that the dependence u d is about right. 
Then I had received some useful advice at my Institute and in one 
of the evenings I carried out a series of measurements, the results 
of which are presented in Fig. 1. 

Following is a short description of the "technology" of the 
experiraent which any reader may carry out by himself at his kitchen. 
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having only a stop-watch and a couple of hours. To measure 
distances I had drawn by ball-point pen a 1 cm grid at the bottom 
of a white enameled kettle 20 cm in diameter. The grid did not 
change in the least during two hours of experiments with water at 
80-90®C. The kettle was in a water bath. The bath was composed 
from a wide frying pan, the kettle standing on three small pieces 
of a wooden rod 1 cm high and the depth of the water (at 90®C) 
was about 2 cm. All the construction was on a slow constant gas 
fire. On a nearby stove was a large teapot with water of about 
the same temperature on a small fire. The teapot water was used 
to add water into the pan (evaporation!) and to change its depth 
in the working kettle. Temperatures were measured by a laboratory 
mercury thermometer. Inside the kettle the mean water temperature 
was 83 + 1®C for all measurements. It was changing little during 
measurements and so if one is interested only in the check of the 
^ ^ d dependence one may not have the thermometer but should keep 
all the fires constant. The depth of the water in the kettle was 
measured by a ruler^ The choice of tracers was a problem at first. 
Meanwhile I found that almost any dry organic powdered material, 
e.g., black pepper, may serve as a tracer, because becoming wet it 
has practically neutral buoyancy. In the experiment described I 
used powdered (by myself) tea and a dry red wild rose berry (also 
powdered) . 


I 
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Most of these particles were sitting at the bottom; part were 
on the surface but some were transported within the water. Their 
path was measured by an eye on the coordinate grid and the time 
of the rectilinear parts was measured by the stop-watch. For 
the results presented in the figure, there were eight layers with 
a depth from 2 to 10 cm. At each layer there were 35 individual 
measurements of path and time. During the time of the measure-*- 
ments at each depth a layer of water about 3 mm thick was evapora- 
ting and the horizontal scale of the points at the figure reflects 
this fact. 



I 


The dispersion of individual velocity measurements around their 
means is in the limits of 10-15% which corresponds, I have said, to 
the accuracy of visual measurements. However, an accuracy of 

each individual measurement seems to be better, therefore these 10-15% 
reflect evidently a natural dispersion of velocities at convection. 
Through all the points with bars a direct line could be drawn. Some 
deviation of the last points might be due to two reasons: (i) an 

increase of heat losses through the kettle wall, and (ii) an under- 
estimate by eye of the vertical component of the path at larger depths 
of the water layer. 

The value of the proportionality coefficient between u and d 
is about 0.2 sec*"^ as follows from the data of Figure 1* For t = 83®C 
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the dynamic viscosity M = 3.5*10“^ P. The values of a and Cp depend 
for water only very weakly on temperature. Knowing u/d = 0.2 sec“^ 
we can estimate the heat flux by the formula f = 12 ycpU^/agd^. it 
is 0.2 W/cm^. The total flux into the kettle is then about 60 W, 
which looks a reasonable value. 

The knowledge of u and f allows one to estimate the basic 
similarity criteria. Let us do this for the shallowest depth of 
2 cm. Then the Reynolds number Re 300, the Peclet number Pe 600 
and the Rayleigh flux number (see eqn. A2*3) Ra^'vlO^. The convection 
is rather irregular though it is not possible to consider the turbul- 
ence as fully developed as, e.g. , in the atmospheric surface layer. 

In fact, knowing the heat flux density and the efficiency y = d/H 10"*^ 
one can estimate the specific rate of the dissipation e = //pH = 10^^ 
cm^sec^^. Using [24] and [33] one finds = n/d ^ 20, which seems too 
small for the existence of a regime of a developed turbulence. 

In concluding this section, it is repeated once. more that the 
exact formulation of the conditions of the validity of the formula 
[1] or [1*J in the sense of the value of the Reynolds number (or the 
Rayleigh number ““ see Appendix 2) is still to be found. Apparently 
new theoretical and extensive experimental studies have to be performed. 
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6. Gravitational convection 

A number of geophysicists believe that the convection in the 
Earth's mantle is caused by differentiation of the mantle's material 
when a heavier fraction is descending and a lighter one is arising 
(see, e.g., Artyushkov, 1968, Sorokhtin, 1974). How the differentia- 
tion is proceeding at a molecular or at some macroscopic level can 
only be speculated upon. It is believed that this process is taking 
place at the mantle- liquid core interface. In such an uncertain 
situation only the simplest phenomenological approach is justified. 
Such an approach will be tried here. 


Equations describing the gravitational or density convection 
may be written in the Boussinesq-like approximation as follows: 


dv _ _ 
dt 


4 - g — + ^ 


[34] 


= 0 , 


[35] 


^ + kj3Ap' 


[36] 


Here is a given field of the mean density, p' is the deviation 
of the density from the mean due to the differentiation causing the 


1 
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convection, is the rate of the density differentiation in the 
volume if such a process is taking place, kj^ is a diffusion coef- 
ficient of the density deficit. It could probably have a sense 
similar to the filtration coefficient of a liquid in a porous media 
(see, e.g., Barenblatt et al. , 1972). 

The energy equation [6] should be added here together with the 
equation of state p=p(T). If the differentiation is taking place 
at some surface (at the layer's lower boundary, say) then it may be 
described by a value r^- with the dimension ML i.e., of the 
density flux. At mathematical modelling it is equivalent to the 
following condition 


rt = kj)9p ' /Sz <at z = 0, say, [37] 

and other boundary conditions are the same as for the thermal 
convection. 


The total flux of the density deficit at a level z consists 


of two parts: 


M (z) 


z r 

rt + 

« 

0 


p ^- ( z)d z . 


[38] 


The flux plays the role analogous to the heat flux in the thermal 
convection. To make the analogy quite clear, note that the heat flux 
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can be connected with the density deficit flux M in the following 
way ; 


-1 


M = p'w' = ap T'w' = acp 


[39] 


where the overbar means some average. The setting of this relation 
allows one to write down at once the formula for the mean velocity 
of the gravitational convection, substituting in [1] M for a//cp: 


u ^ 


Md2 % 




[40] 


pv 


The specification of a numerical coefficient in this formula 
can be done by analysing the energetics of the process as for the 
thermal convection. The flow with the density deficit flux M 
rising to the height d in the gravity field g is releasing the 
potential energy with the rate M gd in the unit column of height 
d (in fact, gp’ is the force on the unit volume, gp'w is the power 
of the force, gp’wd = A?gd is the total power released in the unit 
column. This power is spent on the rate of generation of the 
kinetic energy of the convection which should be equal in average 
to the rate of its dissipation for a steady process. We have 
already estimated in Section 2 the dissipation of kinetic energy 
in the unit column for the plane cell as 


i 


..a. 
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\ Pod-32v^/d2 ^ 


Equating this expression to M gd one gets 



[41] 


Here again is evidenced the analogy with the structure of the 
velocity field for scales less than the viscous microscale. The role 
of the dissipation rate e here plays the value of ^ s/Pq 


Note that from the equality ^ i ^ gd it follows at once that 
the efficiency of the gravitational convection is unity, if M = const (z) 
which is right for the differentiation only at the lower boundary.* 

The substitution into this equality of the density deficit flux for 
the thermal convection from [39] gives i “ ( gd/Cp)/= (d/H)/. It 
looks like a simplest derivation of the formula for the efficiency 
of the thermal convection y = d/H. 


* If the separation of the material takes place only within the layer 
then the total power released at the convection will be equal to 


d 

' ^ 

• Z 1 

Q - 

g M(z)dz = 

gdz 


o 

0 c 


p^(.(z*)dz' 


For g = const and pj- = const Q = % M (d ) gd , where M (d) = p^d. 

The efficiency of such a process is twice less than for the differentia- 
tion at the lower boundary only. A general case may be considered 
similarly to the efficiency of thermal convection (see Hewitt et al., 1975). 
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Therefore, at convection the power l\J gd is transformed at the 
end into heat and ensures the heat flux which could be measured 
at the upper surface. Because in the formation of the observed 
heat flux purely thermal sources can also take part then a constraint 
on the flux of the density deficit follows from here. 

M < //gd . [42] 

This inequality is of importance for the consideration of 
processes in the Earth's interior (see following Section 7). 

The formula for the velocity scale [40] can be also obtained by 
the similarity and dimension arguments. The simplest conditions 
to get it are for = 0, Re << 1 and << ud. The last allows 
us to neglect the diffusive transport of matter- Then in the vorticity 
balance equation obtained by taking curl of [34] we will get a single 
dimensional combination g/vp^ = g/y with the dimension in 

the bouadary condition [37] there is the f lux W of dimension ML“^T”^ 
and there is the vertical scale d. Using these three values one 
constructs the unique velocity scale which is [40]- The time scale 
does not again depend on the geometry scale d: 

T - d/V ^ (y/g 


[43] 
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For the total kinetic energy of the volume d one can get 

E = c [44] 

2v 

where again c ^ 1/32 or 1/72 depending on the two^ or three-dimensionality 
of the convection. It is evident that this formula is similar to the 
corresponding formula [20] for the thermal convection. 


The found analogy allows one to model the gravitational convection 
using the thermal one. Then the role of the Prandtl number Pr = v/lc 
would play the Schmidt number Sc - v/d^ and the role of the Rayleigh 
flux number defined as 





[45] 


would play the number 


Ra = 



[ 46 ] 


A simple example of gravitational convection from everyday exper- 
ience is the motion in a liquid where gaseous bubbles are forming, 
such as in the glass with gasified mineral water or with any other 
bubbling liquid. However, the visualization of these motions and 
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devising any quantitative measureinents proved to be a hard problem 
for me. After many trials and wasting many bottles of the mineral 
water, I finally was able to invent a simple and quick experiment 
based on the check of eqn. [43], i.e., the independence of convection 
time scale on the depth of the liquid. 

As the working fluid I used the mineral water, *‘Moskovskaya” 

(from a drill hole within Moscow City) , which was poured into a 
transparent glass flask of a parallelpiped shape with sizes 95 x 79 
X 37 mm. It had been noted that adding small particles (of powdered 
black pepper*) increase strongly the bubble formation and the in- 
tensity of motions. The marked release of gas takes place in these 
conditions for several hours. The experiment consists of the follow- 
ing procedures. First I put into the flask by a pipelet the black 
pepper in water (about 0.5 cm^). Then the mineral water v;as poured 
into the flask first up to the line of 2 cm. All this was allowed 
to settle for a few minutes to cease the motions caused by pouring. 
Then to the surface of the water a droplet of dye (alcohol solution 
of brilliant green) was introduced by the pipelet. In the water 


* Originally the powdered pepper was thought to be used as a tracer; 
however, after some observations I noticed that the particles’ upward 
motion was caused mainly by gas bubbles sitting on the particles or 
caught by them and the downward motion originated as from an elastic 
strike of the rising particles on the surface film of the water. 
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the droplet formed at once a little cloud from which dye threads or wlsns 
were pulled out. For control such a droplet was introduced into 
an ordinary water where it spread much slower and the dye remained 
mainly in the upper layer of the water. Two typical times were 
measured, the time when a dye thread first touched the flask 
bottom, and the time X7hen many threads spaced more or less 
uniformly were touching the bottom. If the first time may be 
determined more or less distinctly, the second time, T^, is deter- 
mined rather subjectively. The results of two series of such 
measurements are presented in the following table. 


TABLE 


d, cm 1 

Tl, 

sec 

T 

sec 

2 

12 

1 

16 

25 

35 

3 

15 

— 

25 

— - 

4 

— 

11 

— 

30 

5 

1'^ 

— 

25 

— 

6 

— 

15 

— 

30 

7 

12 

16 

30 

30 

8 

— 

14 

— 

25 

9 

15 

15 

30 

25 


The data of the table show that both times may be considered 
as independent on the depth of layer of mineral water. The mean 
velocities of motions calculated from the data of the table change 

from about 1.4 mm/sec to 6 mm/sec and are presented in a lower 
part of Figure 1. 
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For the shallowest layer the Reynolds number is about 30 
(v 0.01 cm^/sec), and for the deepest one it is about 500. Though 
these numbers markedly exceed unity, the experience obtained in 
studying the thermal convection (see Section 5) tells us that here 
some statistical regularities valid for Re << 1 may still be 
operative. Therefore, the absence in our measurements of any 
systematic dependence of the characteristic times of motions on 
the layer depth may be considered as an argument in favor of the 
theory proposed here. 


7. Applications to the mantle convection 

Here a brief discussion will be presented of the developed 
theory to the Earth’s mantle. Because here (following McKenzie et al., 
1974) parameters of the medium are constant, then it may be considered 
only cs the first rough approximation of the convection in the mantle. 

In reality the viscosity depends strongly on the temperature and 
pressure (see Carter, 1976) which exerts important influence on the 
flow pattern. In particular the computations of the thermal convection 
by Houston and DeBremaecker (1975), carried out with the Herring- 
Nabarro viscosity depending exponentially on temperature, give an 
appreciable intensification of the convection in the regions of de- 
creased viscosity and its attenuation there where the viscosity is 
large. Nevertheless these computations and experiments by Booker 
(1976) show that the flow patterns and intensity of the motion are 


1 1 1 I 1 / I 
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not too drastically different from the convection with constant 
parameters of the material. Therefore our formulas are giving right 
orders of convective velocities if one uses some effective vis- 
cosity value. However, this question requires additional studies, 
laboratory and numerical ones. 

The values of thermal convection velocities of order 1 cm/yr 
obtained here (after McKenzie et al. , 1974) seem to be insufficient, 
because many lithospheric plates are moving several times faster. 

If one takes into account that the plates are moving as a whole, 
dragging each other, the oceanic plates are diving under continental 
ones, then one would feel safer if the mantle’s motions would have 
velocities, say, of order 10 cm/yr. 

The structure of the formula [1] shows that this may be reached 
by an increase of the coefficient of thermal expansion a and/or by 
a decrease of the dynamic viscosity y. Hewitt et al. (1975) note 
that the value of a is rather uncertain and could be, in principle, 
increased by an order of magnitude which would increase the velocities 
by a factor of 3. However, this would also mean that the efficiency 
Y ^ d/H 1. But at d E some other complexities arise (see their 
discussion by Hex^itt et al. , 1975) leading to possible inapplicability 
of our simple considerations. One should not exclude that the value 
of viscosity V = 2* 10*“^^ m^sec**^ (or y = 7 . 4* 10^^ P) adopted by 


I I 1 1 I I \ 1 
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McKenzie et al. (1974) and here is also considerably overestimated 
(see also Carter, 1976). So it looks like it is possible to increase 
u up to 10 cm/yr even for the purely thermal convection in the 
upper mantle. 

Let us discuss now the possible gravitational convection. Here 
the whole mantle, including the lower one, will be considered, i.e., 
the layer of some 3000 km thickness, due to the belief that the 
differentiation goes on the mantle-core interface. Not going into 
details of the geophysics of these questions, we present here only 
some estimates of the intensity of such a conveption and point out 
some constraints. 

Consider first inequality [42]: ^ < //gd. Let f = 0.06 W/m^ 
and d = 3000 km = 3*10^ m. Then ^ < 2*10“5 kg/m^sec - 6*10“^ kg/m^yr. 

For the period of time the density of the mantle material will be 
changed by 

We neglect further the non-uniformity of the differentiation 
rate in the process of the Earth's evolution. Due to the models 
of the evolution by Keonjan Monin (1975) and Monin and Keonjan 
(1976) this rate for the last four billion years changes less than 
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twice. Then taking the present value of f and T^ = 4 Aeons, one 
gets from [47] that Ap < 100 kg/m^ =0.1 g/cm^. This value could 
be somewhat increased if one assumes that some part of the heat 
released at gravitational convection goes to the support of the 
differentiation reactions which are, evidently, endothermic. Then 
instead of [42] we should write 

^ gd < / + Q , [48] 

where Q is the heat spent in a unit column for the support of the 
reactions. Nevertheless, the value Ap 1 g cm"*^ adopted by 
Artyushkov (1968) and Sorokhtin (1974) seems to be too high not 
only from the point of view of the constraints [42] or [48], but 
also regarding estimates of the energy released at the gravitational 
differentiation. In fact, due to Monin and Keonjan (1976) and 
several other models the total energy released at the process is 
of order 1.5*10^^ J for the whole Earth’s history. If all this 
energy would be brought up to the surface uniformly, then the geo- 
thermal flux would be of the order G.2 Wm“^, i.e. , thrice the present 

value. The excessive energy can go only to the heating of the 

25 

core. If the core mass is of the order 10 kg then the heating 
for 4 Aeons would be about 2000 K. It would be less if part of 
the energy would go for the support of the differentiation reactions. 
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For Illustrative purposes we present now estimates of mean 
velocities for the gravitational convection. Let M ^ 2*10 
kg/m^sec. For the dynamic viscosity of the (lower) mantle p ^ 10^^ 
kg/m sec = 10^^ P (see McKenzie et al. , 1974) one gets from [41] 
that u 'v 0.1 cm/yr. Monin and Keonjan (1977) assumed the value of 
p by three orders less, as a representative for the whole mantle. 
Then u 3 cm/yr. Not being a specialist on the Earth’s interior 
geophysics, I end the discussion at this point. 



Discussing theories developed here, we have already referred 
to the similarity theory for circulation of planetary atmospheres 
and to Kolmogorov’s theory of turbulence in the viscous range of 
scales. The last is found to be the direct analog of the viscous 
fluid convection studied here. Therefore, we may consider that the 
thermal and gravitational convection at the not too large Reynolds 
numbers and turbulence in the dissipation range are forming a family 
of forced viscosity dominated flows. 

However , there are more common features among all the afore- 
mentioned kinds of flow. The first is that the total kinetic/ 
energy of a general circulation and of convection does not depend 
on the total mass of the atmosphere or of the convec ting fluid. 





I 




For convection this property is evident from eqn. [20] for the 
thermal case and eqn. [44] for the gravitational case. The kinetic 
energy of the circulation of a planetary atmosphere is equal (up to 
a multiplier or a non-dimensional function of angular velocity) to 

E , [49] 


where cf is the Stefan-Boltzmann constant, <1 = cjf^(l-A)/4 is the mean 
rate of solar energy reaching the planet accounting for its albedo A, 
r is the planetary radius (see Golitsyn, 1970, 1973). Independence 
of the kinetic energies on the mass of flows is a more general feature 
which might be used for a classification. 


However, many types of flow do not have this property, but 
nevertheless, they have another more general feature which we shall 
discuss now, starting from circulations. The eqn. [49] after some 
simple transformations can be written as 



[49’] 


where K = Cp/c^, Q = 4irr2q - the total power of the solar energy 
assimilated by a planet, /= t( is the sound velocity, 

Te - (q/o^)^ is the temperature of the eqiulibriura radiation of 
the planet. The quantity = r/c^ is the time for propagation of 
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a perturbation in the atmosphere in global scale. As is knoxm 
(Landau and Lifshitz, 1959, §48) the sound velocity is character- 
istic for reaching the local thermodynamic equilibrium. Therefore 
up to a multiplier of order unity the total kinetic energy of 
circulation is 


E Q t49**] 

i.e., it is equal to the total power of the radiation assimilated by 
a planet times the time of perturbation relaxation in the global scale. 

Expression [44] for the total kinetic energy of the gravitational 
convection E % A/ gd^/v has exactly the same structure. The combina- 
tion q g = gd is the power of this convection in a unit column. The 
total power of the convection in the whole volume (or the rate of 
potential energy release in the volume) is Q = q d^ = Af dg^. From 

o 8 
9 

the other hand, d /v = is the viscous relaxation time in the 
volume with characteristic size d. As a result 

Eg ^ gd^/v = 144'] 

The peculiarity of the thermal convection is that only some 
part of the heat power brought . in a fluid is spent; on the generation 
of kinetic energy of convective motions. The part determined by 
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the efficiency coefficient of the convection is equal to y - agd/cp 
= d/H. Accounting for this circumstance, eqn. [20] can be rewritten 
as 

-i ^ . (20’1 

where F = fd^ is the total heat flux introduced into the fluid. The 
same structure is having eqn. [27] for the kinetic energy of a volume 
of turbulent fluid with size d < n . 

The similar form may be given also to the expression for the 
kinetic energy of a volume of locally isotropic and homogeneous 
turbulent flow relative to its closest volumes with the size of 
n < d < Lg, where is the turbulence external scale, i*e,, for 
d in the inertial interval (Kolmogorov, 1941). The energy of such 
a volume with the mass M = pd^ is 

E = Qd/(Ed)^^^ , [50] 

where Q = Me is the total energy rate brought into the volume. 
Because in the inertial range of turbulence (ed)^^^ V the 
relative rms velocity for tx^o points separated by the distance d 

E Q-d/v , [50'] 


then 
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/ 

where is the characteristic lifetime of eddies with scale d. 

Note, hovzever, that the locally isotropic and homogeneous 
flow in the inertial range does not belong to the class of flows 
whose hinetic energy (in our relative sense) does not depend on 

the fluid mass. Because e = Q/M the relative kinetic energy, as 

— 1 / * 
follows from [50], is proportional to M ' • Analogous dependence 

of the total kinetic energy from mass has been obtained in some 

models of general circulation, considered as a large scale convection 

on non-rotating planets (Gierasch et al. , 1970, Leovy and Pollack, 

1973, Burangulov and Zilitinkevich, 1976). In these models one 

may also obtain formulas of the type of [50]. In Appendix 3 we 

consider a turbulent convection of a fluid in a bounded region. 

This kind of convection also belongs to this type of flow. 

Together with the aforementioned models of circulation and turbulent 

flow in the inertial range it forms a family of flows with mean 

velocities proportional to (sd)^^^. This group of flows may be 

called a family of forced turbulence dominated flows or inertial, j, 

because the inertial non-linearity and resulting turbulent mixing 
determine their structure and intensity. 

We see that for quite a number of forced geophysical flows • 

their total kinetic energy is determined by the product of total 
power brought into the fluid and of characteristic relaxation time. 


f 


1 


I 


I 
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Note that in all the cases considered here this time is always 
the smallest from all the times which can be constructed from ex- 
ternal parameters in the equations. It is true, however, that this 
suiallest time is usually only one if we believe in the validity of 
corresponding self-similarity hypotheses neglecting various external 
parameters* This allows one to propose the following approximate 
rule which could be called a ’’principle of the fastest response”. 

The ’’principle” says that the kinetic energy of a forced steady 
flow of a fluid system is of the order of total power brought into 
the system times the shortest relaxation time characteristic to the 
system. 

If one is not using the similarity theory then this ’’principle” 
allows one to write at once the expression for the total kinetic 
energy of motions. It was this '’principle” noted by me for the general 
circulation x^hich has been used in 1970 (Golitsyn, unpublished) for 
obtaining an expression of the kind of eqn. [20], but without account*- 
ing for the convection efficiency y (it looks like nobody at the time 
had any idea of it). I was getting then V ^ 30 cm/yr for the mantle 
and the excess was mainly not due to not accounting for efficiency 
but due to not realizing that the numerical constant c in [20] should 
be about 1/30 (or even less) . 

We see that obedience of the flows considered here to the 
’’principle of the fastest response” is their most general property* 


I i I i I I k - 

-56- 

However, to drive away an Impression on the universality of the 
principle I want to present an example of a system where it gives, 
at bestj an estimate from below of the kinetic energy. This is 
the circulation in atmospheres of large and fast rotating planets 
Jupiter and Saturn. A detailed discussion of their circulation 
was given by Golitsyn (1970, 1973). It looks like the fast rotation 
is a factor strongly stabilizing large-scale motions and not allow- 
ing the system to relax by the fastest way. 

Results of this section are mainly of methodical character. 
However, it seems that the ^'principle of the fastest response*’ may 
also have a heuristical value as it has had for this paper for 
which it served as a first impulse. 
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APPENDIX 1 

SOME GENERAL FORMULAS TAKING INTO ACCOUNT THE EFFICIENCY OF CONVECTION 


The concept of a convection efficiency measured by the 
coefficient Y = d/H = agd/cp, suggested by Lliboutry (1972) and by 
Hewitt et al, (1975) seems to be an important achievement of hydro- 
dynamics the use of which has not yet spread wide enough. In this 
connection it should be noted that the quantity y implicit 

form enters in fact many important formulas and definitions of hydro- 
dynamics. )Io help clarify this concept I wish to rearrange several 
formulas into a form where the efficiency y will be present explicitly. 

First of all two formulas where there is a temperature difference 
AT. Introduce the entalpy difference Ae = CpAT. Then the Rayleigh 
number can be transformed as follows : 

Ra = = i Y lie [A.1.1] 

KV H KV KV 

Similarly, the thermal Rossby number is transformed: 

Ro„ = = Y [A. 1. 2] 

x^here 0 is the angular velocity of rotation of a vessel with a 
lateral scale b * Convection there from the point of view close 
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to the one used here has been consider by the author (1977b). 

Finally, the non-dimensional height in the Monin-Oboukhov 
similarity theory of the stratified turbulent boundary layer 
(see Monin and Yaglora, 1965, §7). 

? = Z = Y(z) J2L [A.1.3] 

pcpu3 pu3 

where y(z) = ^/H and is the friction velocity. Though y does 

* 

not play here the role of the efficiency, its presence in these for- 


mulas seems to be instructive. 
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APPENDIX 2 

RELATION OF THE DEVELOPED ELEMENTARY THEORY OF CONVECTION WITH USUAL 

REPRESENTATIONS 


Usual similarity criteria for the thermal convection in the 
existing literature are the Rayleigh and Prandtl numbers. The 
first of these have not yet been present in the theory here. As 
we will see later, this happened because we were not considering 
the structure of thermal boundary layers. Let us show how the 
Rayleigh number is emerging here and how it is connected with the 
similarity criteria which have already been introduced. 


In the non-dimensional system (15*)- (17’) there are three 
similarity criteria: (i) the Reynolds number 



T 

V 





[A. 2.1] 


(ii) the Peclet number 



[A.2.2] 


and (iii) the efficiency coefficient y At the same time from the 
dimensional parameters present in [15]-[17] one may construct the 
Rayleigh number defined by the heat flux: 
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. ssOH , lil . . 

j pc k^v |ik^ k^pH 


Ra 


[A. 2.3] 


Comparing this expression with [A. 2.1] and [A. 2.2] we get 


k , 


Re = Ra^^/Pr , 


[A.2.4] 


Pe = 


[A. 2. 5] 


To connect the Rayleigh flux number with an ordinary Rayleigh 
number [A. 1*1] one has to know the dependence of the Nusselt number 
(see eqn. [30]) on the Rayleigh number. The study of the dependence 
Nu(Ra) is the usual goal in the laboratory convection experiments. 

It may be also obtained from theoretical considerations of the energy 
and vorticity balance in the boundary layer (see Turcotte and Oxburgh, 
1967, McKenzie et al. , 1974). The measurements and the considerations 
produce the dependence 


Nu ^ 0 Ra 


1/3 


[A.2.6] 


where g is some numerical coefficient (about 0. 1-0. 2 in dependence 
on” the problem conditions) . Using [A. 2-3], [A. 1.1] and [A.2.6] one 

gets 


Ra^f = gRa^^^, or Ra = (Ra^/$)^/^ . 


IA.2.7] 
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Comparing [1'] with [A. 2.3] or [A. 2. 7] we may obtain 

I = I 

Expressions of this kind but without the numerical coefficient 
have been obtained by Turcotte and Oxburgh (1967) and McKenzie et al, 
(1974). Using the Rayleigh flux number one can express also the 
temperature difference in the main body of the fluid: 


6T 


. ^ <32 


% 

Oigd^ 


[A.2.9] 


i 


Related to the total temperature difference AT causing the 
convection the difference 6T is equal to 


AT 


(32 Raf) 


Ra 


( 326 )^ 


Ra 


-1/3 


% 


32”^3 3/4 


Ra 




■/ 


[A. 2 .10] 


Therefore with the increase of the Rayleigh number the temperature 
profile would appear more and more like isothermal one in the bulk of 

fluid comparing with temperature drops in the thermal boundary layers. 
Laboratory experiments and numerical computations (see, e. g. , McKenzie 
et al., 1974) show that for a developed convection the main temperature 
changes are in the boundary layers of about %AT in each (upper and 


lower ones). Let us define the thickness of the boundary layers 
6 using the boundary condition [17] in such a way that 


26 

Comparing this expression with [30] we get 


[A.2.11] 


Nu = d/26 , 


from where accounting for [A. 2. 6], it follows that 


6 (%-e)d Ra d Ra^^"^ 


[A. 2.12] 


If at the lower boundary we know the heat flux but the total 
temperature change AT is not known, we may estimate it from [A. 2. 9] 
taking into account [A. 2, 12] as 



[A. 2.13] 


The expressions of the type of [A. 2.12] and [A. 2.13] are al'^o 
in agreement in McKenzie et al. (1974) and results of their com- 
putations agree well with the dependencies 6 a; f and AT ^ 
following from the expressions. 
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The last thing which is useful to estimate is the ratio of 
temperature gradients in the main interior of the fluid — and 
in the boundary layers — • follows from [A, 2. 10], [A.2\12] 

and [A. 2. 13] that 


V^T 


% 


32 % 

33/A 



6 

d 



[A. 2. 14] 
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APPENDIX 3 

TURBULENT CONVECTION 


At last we consider the turbulent convection when the Reynolds 
number is very large and the direct role of viscosity in the bulk 
of the fluid is not important. Exact conditions for the existence of 
such a regime in terms of external parameters or similarity criteria 
have yet to be established mainly experimentally, in the laboratory 
or by computations. Though this regime comes off the frame of the 
basic subject of this paper, nevertheless its consideration is justified 
here from the point of view of Section 8, because the regime as we 
will see later is the closest analog of the regime of developed turbul^ 
ence, described by the Kolmogorov-Oboukhov theory. This adds to our 
scheme of classification of flows by one more object. 


Based on the vorticity equation [15] we see that for a steady 
case the main balance is between inertial and buoyancy terms. 
Estimating the entalpy gradient from the energy equation [16] and 
substituting it into [15] we obtain the following estimates of scales 
of velocity and entalpy gradient (again at the main body of fluid). 



From the last formula an expression follows for the temperature 
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difference In the bulk of fluid with height of order d: 


1 / 3 , 


\2/3 
P ^ 


(0 


[A. 3. 3] 


The expressions [A. 3.1] and [A. 3. 3] will be identical to the 
expressions for variances of velocity and temperature for the free 
convection in the atmosphere if one substitutes there the running 
height z for the layer depth d and remember that for a gas a = 1/T, 
where T is a characteristic temperature. Let us remember that the 
expression for the velocity variance in atmospheric convection has 
first been obtained by Prandtl (1932), and the expression for the 
temperature variance has been obtained by Oboukhov (1960), who being 
unaware of the Prandtl paper has also obtained an expression of the 
[A. 3.1] type and’ several other useful expressions. Both these scientists 
were using similarity and dimensional arguments and here the scaling 
analysis is used. 

1/3 

Because here eqn. [A. 3.1] could be written as u (^d) . This 

clarifies the analogy between regimes of turbulent convection and of 
locally isotropic and homogenous turbulence. We see that Prandtl 
(1932) , nine years before Kolmogorov (1941) was not too far from dis- 
covering the laws of local structure of turbulence. 


If we use the scales u and Ve for non-dimensionalizing the 
equation system [15]- [17], then we will have as the similarity 
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criteria the Reynolds and Pdclet numbers 

■ T “ f ■ [A. 3. 51 

As another similarity criterion for convective flows we may again 
take the Rayleigh flux number defined by [A. 2. 3]. Here its relations 
with the other criteria differ from the previous case and are of the 
form 


Re = ^ Pr“2/3 

Pe = Ra^^^ Pr^^^ 

/ 


[A. 3. 6] 
[A. 3. 7] 


The temperw^ture difference in the bulk of fluid is in these terms 
equal to 


6T ^ Ra 


•/ 


2/3 (k2y)2/3 


agd' 


[A. 3. 8] 


and the velocity scale 


u ^ (Ra^ Pr)^/3 I 


[A. 3. 9] 
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* 

\ 

The structure of the thermal boundary layer remains the same 
as in the viscous case considered in Appendix 2 because the dependence 
Nu ^ Ra is valid here too. The last formula here is the 

ratio of temperature gradient in the bulk of fluid and in the boundary 
layer: 

ViT -1/3 

^ ^ (Ra^ Pr) [A.3.10] 
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Figure 1 


The dependence of rms velocity u of the fluid on its depth. 

Upper solid dots with standard deviations (bars) are for the 

thermal convection. The lower part is for the gravitational 

convection. d/T-, is the maximal velocity, u = d/T is 

1 1 m m 

the mean velocity, dots and crosses are for two runs of the 
experiment according to Table 1 (see also text). 




